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Abstract 

• We improve the theoretical estimates of the critical exponents for the three-dimensional XY 

O 

universality class, which apply to the superfluid transition in 4 He along the A-line of its phase 
diagram. We obtain the estimates a = -0.0151(3), v = 0.6717(1), n = 0.0381(2), 7 = 1.3178(2), 
(3 = 0.3486(1), and 5 = 4.780(1). Our results are obtained by finite-size scaling analyses of high- 
statistics Monte Carlo simulations up to lattice size L = 128 and resummations of 22nd-order 
high-temperature expansions of two improved models with suppressed leading scaling corrections. 
We note that our result for the specific-heat exponent a disagrees with the most recent experimental 

o ; 

estimate a = —0.0127(3) at the superfluid transition of 4 He in microgravity environment. 
PACS numbers: 05.70.Jk, 64.60.Fr, 67.40.-w, 67.40.-Kh 
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I. INTRODUCTION AND SUMMARY 



The renormalization-group (RG) theory of critical phenomena classifies continuous phase 
transitions into universality classes, which are determined only by a few global properties 
of the system, such as the space dimensionality, the nature and the symmetry of the order 
parameter, the symmetry-breaking pattern, and the range of the interactions. Within a given 
universality class, the critical exponents and scaling functions describing the asymptotic 
critical behavior are identical for all systems. The three-dimensional XY universality class is 
characterized by a complex order parameter and symmetry breaking 0(2) = Z 2 <g)U(l) — > Z 2 . 
An interesting representative of this universality class is the superfluid transition of 4 He along 
the A-line T\(P), which provides an exceptional opportunity for a very accurate experimental 
test of the RG predictions, because of the weakness of the singularity in the compressibility 
of the fluid and of the purity of the samples. Exploiting also the possibility of performing 
experiments in a microgravity environment, the specific heat of liquid helium was measured 
up to a few nK from the A-transition. 1 The resulting estimate of the specific-heat exponent, 
obtained after some re-analyses of the experimental data, 1 ' 2 ' 3 is 3 



Other experimental results at the superfluid transition of 4 He, and for other physical systems 
in the XY universality class, are reported in Ref. 4. 

On the theoretical side the XY universality class has been studied by various approaches, 
such as field-theoretical (FT) methods and lattice techniques based on Monte Carlo (MC) 
simulations or high-temperature (HT) expansions. A review of results can be found in 
Ref. 4. Accurate estimates of the critical exponents were obtained in Ref. 5 by combining 
MC simulations and HT expansions of improved Hamiltonians. This synergy of lattice 
techniques provided the estimate 



which is substantially consistent with the experimental result (1), although slightly smaller. 
Other results obtained from FT calculations, 6 ' 7,8 MC simulations, 5 ' 9 and HT expansions 10 ' 11 
are less precise and in agreement with both estimates (1) and (2). 

In this paper we significantly improve the theoretical estimates. This allows us to check 
whether the small difference between the theoretical and experimental estimates (1) and 



a = -0.0127(3). 



(1) 



a = -0.0146(8), 



(2) 
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FIG. 1: Summary of our results for the specific-heat exponent a. Shorthands are explained in the 
text. The coloured region corresponds to our final estimate a = —0.0151(3). 

(2) disappears after a more accurate theoretical analysis, as recently suggested in Ref. 12. 
We again follow the strategy of Ref. 5, considering two classes of lattice Hamiltonians, 
the 4 lattice model and the dynamically diluted XY (ddXY) model. They depend on an 
irrelevant parameter, A and D respectively, which can be tuned to suppress the leading 
scaling corrections, giving rise to improved Hamiltonians. We improve the finite-size scaling 
(FSS) analysis of these models by significantly increasing the statistics (by approximately a 
factor 10) and simulating larger lattices (up to lattice sizes L = 128). The precision of the 
data allows us to observe the expected next-to-leading scaling corrections, and therefore to 
have a much better control of the systematic errors. Moreover, we extend the HT expansion 
of the susceptibility and of the correlation length in the 4 and ddXY models to 22nd order, 
i.e., we add two terms to the HT series computed and analyzed in Refs. 5,11. Using this 
bulk of new data and calculations, we performed several analyses, also combining information 
obtained from MC simulations and IHT analyses (IHT denotes the HT expansion specialized 
to improved models). 

In Fig. 1 we show a summary of our results for the specific-heat critical exponent a, as 
obtained from the hyperscaling relation a = 2 — 3z/. We show three sets of new estimates for 
the </> 4 and ddXY models, obtained by different methods. We first report the results of FSS 
analyses of MC data up to L = 128 (FSS). The three reported results for each model are 
obtained from the analysis (left) of a combination of two quantities (temperature derivatives 
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of RG invariant quantities) that does not have leading scaling corrections, (center) of the 
energy density, and (right) of the fourth-order (Binder) cumulant of the magnetization. 
Next, we report the results obtained from the analyses of the 22nd-order IHT series at the 
optimal values of the irrelevant parameters A* and D*, biased using the MC estimates of 
j3 c (MC+IHT). Finally, we show results obtained by requiring the FSS and IHT analyses 
to give consistent results for the exponent v (FSS+IHT). These results come from different 
analyses of the available MC data and HT calculations. Although they are not completely 
independent, their comparison represents a highly nontrivial cross-check of the results and 
of their errors. In Fig. 1 we also show our earlier MC+IHT result obtained in Ref. 5. 

The results shown in Fig. 1 provide a rather accurate estimate of a, which we summarize 
by taking 

a = -0.0151(3) (3) 

as our final estimate. This result significantly improves our earlier estimate (2), but it 
does not reduce the difference from the experimental result (1). Moreover, now the errors 
are so small that their difference appears significant. According to our analyses, values of 
the specific-heat exponent a > —0.014 appear to be highly unlikely. We think that this 
discrepancy calls for further investigations. We mention that a proposal of a new space 
experiment has been presented in Ref. 13. 

We also anticipate our best estimates of the other critical exponents: 



v = 0.6717(1), (4) 

7^ = 0.0381(2), (5) 

7 = 1.3178(2), (6) 

5 = 4.780(1), (7) 

= 0.3486(1). (8) 



Moreover, we obtained an accurate estimate of the exponent u associated with the leading 
scaling corrections, i.e., oo = 0.785(20), and of the exponent uj 2 associated with the next-to- 
leading scaling corrections, i.e., u; 2 = 1.8(2). 

The paper is organized as follows. In Sec. II we define the 4 and ddXY lattice mod- 
els. Sec. Ill is dedicated to a summary of the basic RG ideas concerning FSS in critical 
phenomena. In Sec. IV we present our MC simulations and the FSS analyses of the data. 
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The computation and analysis of the HT expansion are presented in Sec. V. We report 
an analysis of the IHT expansions of improved 4 and ddXY models biased using the MC 
estimates of (3 C (MC+IHT), and a combined analysis requiring the consistency of the IHT 
and FSS results (FSS+IHT). 



II. LATTICE MODELS 



As in Ref. 5, we consider two classes of models defined on a simple cubic lattice and 
depending on an irrelevant parameter. They are the 4 lattice model and the dynamically 
diluted XY (ddXY) model. The Hamiltonian of the 4 lattice model is given by 

= + £ [$l + Kil - !) 2 ] , (9) 

(xy) x 

where (f> x = ((p x ,<f>x ) is a two-component real variable. The ddXY model is defined by the 
Hamiltonian 

n dd = -pY,f*'$v- D 12% ( 10 ) 

(xy) x 

with the local measure 



and the partition function 



(11) 



e- Hdd . (12) 



The parameters A in and D in Hdd can be tuned to obtain improved Hamiltonians 
characterized by the fact that the leading correction to scaling is absent in the Wegner 
expansion of any observable near the critical point. Considering, for instance, the magnetic 
susceptibility Xi the corresponding Wegner expansion is generally given by 

x = Cr 1 (1 + oo,it + a 0<2 t 2 + ... + a 1A t A + a lj2 t 2A + ... 

+h A t 1+A + M 1+2A + • • • + «2,i* A2 + ■■■)> (13) 

where t = 1 — (3/f3 c is the reduced temperature. We have neglected additional terms due 
to other irrelevant operators and terms due to the analytic background present in the free 
energy. 14 ' 15 ' 16 The leading exponent 7 and the correction-to-scaling exponents A,A 2 ,..., 
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are universal, while the amplitudes C, a^j, b^j are nonuniversal. For three-dimensional 
XY systems, the value of the leading correction-to-scaling exponent is A 0.53, 5,6 ' 9 and 
the value of the subleading nonanalytic exponent is A 2 ~ 1.2. 17 In the case of improved 
Hamiltonians the leading correction to scaling vanishes, i.e., a^i = in Eq. (13) (actually 
ai ; j = for all i), in the Wegner expansion of any thermodynamic quantity. 

Improved Hamiltonians belonging to the XY universality class have been discussed in 
Refs. 5,9,11. We mention that improved Hamiltonians have also been considered for other 
universality classes, such as the Ising 18 ' 19 ' 20 ' 21 and Heisenberg universality classes. 22 ' 23 

III. FINITE-SIZE SCALING 

In this section we summarize some basic results concerning FSS in critical phenomena. 
The starting point of FSS is the scaling behavior of the singular part of the free energy 
density of a sample of linear size L (see, e.g., Refs. 15,24): 

.FsingK u h , {M, L) = L~ d F sing (L y *u t: L^u h , {L*Ui}), (14) 

where u t = u±, = u 2 , {u,j} with % > 3 are the scaling fields (which are analytic functions 
of the Hamiltonian parameters) associated respectively with the reduced temperature t 
(ut ~ t), magnetic field H {un ~ H), and the other irrelevant perturbations with < 0. The 
scaling behavior of the interesting thermodynamic quantities can be obtained by performing 
the appropriate derivatives of Eq. (14), with respect to t and H. Since u t and Uh are 
assumed to be the only relevant scaling fields, one may expand with respect to the arguments 
corresponding to the irrelevant scaling fields. This provides the leading scaling behavior and 
the power-law scaling corrections. 

The RG exponents of the relevant scaling fields u t and Uh are related to the standard 
exponents v and 77, i.e., y t = \jv and yu — (d + 2 — rj)/2. The RG exponent 2/3 = — uj of 
the leading irrelevant scaling field w 3 has been estimated by the analysis of high-order FT 
perturbative expansions, 6 obtaining uo = 0.802(18) (e expansion) and uo = 0.789(11) (d — 3 
expansion). Results from lattice techniques are in substantial agreement, see Ref. 4. As 
we shall see later, our FSS analysis provides the estimate u = 0.785(20). Concerning the 
next-to-leading scaling corrections, we mention the FT results of Ref. 17: y± = —1.77(7) and 
t/5 = —1.79(7) (1/421 and y 422 in their notation). Note that, at present, there is no independent 
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check of these results. As we shall see, our new MC and HT analyses confirm that the next- 
to- leading scaling corrections are characterized by a RG exponent y± = —1.8(2). There are 
also corrections due to the violation of rotational invariance by the lattice; the corresponding 
RG dimension is ye = —2. 02(1). 25 We mention that for some quantities, such as the Binder 
cumulant and the ratio = £/L, there are also corrections due to the analytic background 
of the free energy. This should lead to corrections with y 7 = —(2—rj) = —1.9619(2) (obtained 
by using the estimate of rj of this work). Finally, in the case of we expect also 0(L~ 2 ) 
corrections, 26 related to the particular definition (21) of £. 

For vanishing external field H, the behavior of a phenomenological coupling R, i.e., of a 
quantity that is invariant under RG transformations in the critical limit, can be written in 
the FSS limit as 

R(L, (3, A) = r (u t L yt ) + ]T r 3)fc (u t L») u\L^ + ]T r^V*) Ui L* + ..., (15) 

k i>4 

where we have singled out the corrections due to the leading irrelevant operator. The 
functions r (z), r^z), r 3 ^(z) are smooth and finite for z — > and, by definition, u t ([3, A) ~ 
t = 1 — f3/f3 c . In general, all others scaling fields u,i are finite for t — 0. Improved models 
are characterized by the additional condition that u 3 (t = 0) = 0: in this case, all corrections 
proportional to L kys = L~ ku} vanish at the critical point t — 0. In the limit t — > and 
u t L yt ~ tL x l v -> 0, we can further expand Eq. (15), obtaining 

R(L, (3, A) = R* + ct(P, A) tL yt + Ci(f3, A) L Vl + 0[t 2 L 2y \ L 2y \tL yt+y % (16) 

i 

where R* = r (0). As we already mentioned, in improved models all corrections proportional 
to L ky ' A vanish for t — 0. 

Instead of computing the various quantities at fixed Hamiltonian parameters, one may 
study the FSS keeping a phenomenological coupling R fixed at a given value Rf. This means 
that, for each L, one considers @f(L) such that 

R(L,(3 = (3 f (L)) = R f . (17) 

All interesting thermodynamic quantities are then measured at (3 — flf(L). The pseudocrit- 
ical temperature f3f(L) converges to (3 C as L — > oo. The value Rf can be specified at will, 
as long as Rf is taken between the high- and low-temperature fixed-point values of R. The 
choice Rf = R* (where R* is the critical-point value) improves the convergence of/3/ to (3 C for 
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L -> oo; indeed 5 ' 20 /3/ - (3 C = 0(L-V") for generic values of i?/, while /3 7 - & = 0(L _1 /'" _a; ) 
for i?/ = i?*. This method has several advantages. First, no precise knowledge of (3 C is 
needed. Secondly, for some observables, the statistical error at fixed Rf is smaller than that 
at fixed f3 = f3 c . 

Typically, the thermal RG exponent y t — \/v is computed from the FSS of the derivative 
of a phenomenological coupling R with respect to (5 at f3 c . Using Eq. (15) one obtains 



OR 



1 



(18) 



^(0) L yt + r'M Ui(Pc) L yi+y * + n(0) u'm L Vi + .. 

i=3 i=3 

An analogous expansion holds for Sr computed at a fixed phenomenological coupling. The 
leading corrections scale with L y:i . However, in improved models in which u 3 (/3 c ) = 0, the 
leading correction is of order L Vi . Note that corrections proportional to L m ~ yt ps L -2 3 are 
still present even if the model is improved. If one computes the derivative at (3 C , one should 
also take into account the uncertainty on f3 c in the error estimate of v\ therefore, also in this 
case it is more convenient to evaluate Sr at f3f defined in Eq. (17). 

Before concluding the section, we would like to recall three basic assumptions of FSS when 
considering boundary conditions consistent with translation invariance, such as periodic 
boundary conditions: 15 

(a) l/L is an exact scaling field with no corrections proportional to l/L 2 , l/L 3 , etc.; 

(b) the nonlinear scaling fields have coefficients that are L independent; 

(c) the analytic background present in the free energy depends on L through exponentially 
small terms. 

The theoretical evidence for these three hypotheses is discussed in Ref. 15. Under these 
assumptions, there are no analytic l/L corrections. These assumptions can be verified 
analytically in the two-dimensional Ising model, see, e.g., Refs. 16,27,28 and references 
therein. As far as we know, all numerical results reported in the literature are in full 
agreement with assumptions (a), (b), and (c). In particular, there is no evidence of l/L 
corrections to FSS. As we shall see, the FSS analysis that we shall present in this paper will 
provide further support to the FSS assumptions, and in particular to the absence of l/L 
analytic corrections. 
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IV. MONTE CARLO SIMULATIONS AND FINITE-SIZE SCALING ANALYSES 



In this section we present MC simulations that significantly extend those of Ref. 5. The 
statistics are much larger and we consider larger lattice sizes (the largest lattice has L = 128). 
As we shall see, the new MC data allow us to perform a more accurate FSS analysis, 
achieving a much better control of the next-to-leading scaling corrections, and therefore of 
the systematic errors related to the subleading scaling corrections. 

A. Monte Carlo simulations 

We simulated the 4 and ddXY models on simple cubic lattices of size L 3 with periodic 
boundary conditions, at several values of the Hamiltonian parameters A and D, close to 
the optimal values A* and D*, at which leading scaling corrections vanish. Most of the 
simulations correspond to A = 2.07 and D = 1.02, that represent the best estimates of A* 
and D* of Ref. 5. 

The basic algorithm is the same as in our previous numerical study reported in Ref. 5. 
We use a combination of local and cluster 29 updates. We perform wall-cluster updates: 19 we 
flip all clusters that intersect a plane of the lattice. A cluster update changes only the angle 
of the variables. An ergodic algorithm is achieved by adding local updates, 30 which can also 
change the length of the spin variables. We use Metropolis and overrelaxation algorithms 
as local updates, which we alternate with the cluster updates. 

The main difference with respect to our previous numerical work 5 concerns the random 
number generator. Since we planned to increase the statistics by approximately one order 
of magnititude, we decided to use a higher- quality random number generator, such as those 
proposed in Ref. 31. In our MC simulations we used the ranlux random number generator 
with luxury level 2. Its main drawback is that it requires much more CPU time than the 
G05CAF generator of the NAG-library 32 which we used in Ref. 5. In order to get a good 
performance, despite the use of the expensive random number generator, we used demonized 
versions 33 ' 34 of the update algorithms, which allowed us to save many random numbers. 

Most simulations of the 4 and ddXY models were performed at the estimates of A* and 
D* obtained in Ref. 5, i.e., A = 2.07 and D = 1.02. In addition, we perfomed simulations 
at A = 1.9, 2.1, 2.2, 2.3 for the 4 model, and at D = 0.9, 1.2 for the ddXY model. We shall 
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also present some MC simulations of the standard XY model. In total, the MC simulations 
took approximately 20 years of CPU-time on a single 2.0 GHz Opteron processor. 



B. Definitions of the measured quantities 



The energy density is defined as 



1 - - 

E = — fix ■ fiy ■ 



(19) 



The magnetic susceptibility x an d the correlation length £ are defined as 

X 

and 



X/F 



4sin 2 7r/Z/ 



where 



1 /I flnxi 



(20) 



(21) 



(22) 



V W ^ r \ L 

X 

is the Fourier transform of the correlation function at the lowest non-zero momentum. 

We also consider several so-called phenomenological couplings, i.e., quantities that, in the 
critical limit, are invariant under RG transformations. We consider the Binder parameter 
U± and its sixth-order generalization £7 6 , defined as 

<(m 2 ) J > 



U 2j = 



(23) 



where rh = y J2 X fix ^ s the magnetization of the system. We also consider the ratio Rz = 
Z a /Z p of the partition function Z a of a system with antiperiodic boundary conditions in one 
of the three directions and the partition function Z p of a system with periodic boundary 
conditions in all directions. Antiperiodic boundary conditions in the first direction are 

— # — # 

obtained by changing the sign of the term <p x ■ <p y of the Hamiltonian for links (xy) that 
connect the boundaries, i.e., for x = (L,x 2 ,xs) and y = (l,x 2 ,xs). Finally, we define the 
helicity modulus T. For this purpose we introduce a twisted term in the Hamiltonian. More 
precisely, we consider the nearest-neighbor sites (x, y) with x\ — L, y ± — 1, x 2 = y 2 and 
x 3 = 1/3; an d replace the term <p x ■ <j) y in the Hamiltonian with 



4 • Rj y = fi x 1] M 1} cos y? + sin ip) + (0 cos <p - 0« sin <p) , 

10 



(24) 



where R v is a rotation by an angle ip. The helicity modulus is defined by 

1 d 2 \nZ{p) 



T EE — 



(25) 



L dip 2 

The quantities U 2 j, Rz = Z a /Z p , = £/L, and Ry = TL are invariant under RG transfor- 
mations in the critical limit. Thus, they can be considered as phenomenological couplings. In 
the following we will generally refer to them by using the symbol R. Finally, we also consider 
the derivative of the phenomenological couplings with respect to the inverse temperature, 
i.e., 

S* = w (26) 
which allows us to determine the critical exponent v through Eq. (18). 



C. Determination of the critical temperature 

As a first step of the analysis we determine the inverse transition temperature f3 c for 
various values of A in the 4 model and of D in the ddXY model. For this purpose we 
employ the standard Binder-crossing method, i.e., j3 c is determined by requiring 

R([3 C ,L)=R* (27) 

independently of L. Here R is a phenomenological coupling, R* is its fixed-point value, and 
corrections to scaling are ignored. In practice we compute R(j3, L) as Taylor series up to 
third order around the simulation value [3 vun . We choose as [3 rnn our previous best estimate 
of /3 C , see Ref. 5. In the analysis we consider four different phenomenological couplings: 
Rz = Z a /Z p , R^ ee £/L, U4, and Uq. We do not use the helicity modulus because it has a 
poor statistical accuracy for large lattices, which are important to determine f3 c . 

We first analyze the data at A = 2.07 for the 4 model and D = 1.02 for the ddXY model, 
for which we have data with higher statistics especially for the largest lattices. We perform 
fits with ansatz (27) for the two models separately. The two models provide consistent 
results for R*, as required by universality. In order to improve the statistical accuracy, we 
also perform joint fits of the results for both models, imposing the same value of R*. We 
perform fits with various L min and L max (respectively the smallest and largest lattice size 
taken into account), to check stability. In Table I we report our final results, which are taken 
from joint fits of the results for the two models (A = 2.07 and D = 1.02) with 48 < L < 128. 
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TABLE I: Estimates of (3 C at A = 2.07 and D = 1.02 and of the fixed-point value R* for several 
dimensionless quantities. Results of fits of the MC data in the range 48 < L < 128. For details, 
see the text. 



quantity (3 C at A = 2.07 (3 C at D = 1.02 R* 

R z 0.5093835(2) [3] 0.5637963(2) [2] 0.3203(1)[3] 

R^ 0.5093836(2) [3] 0.5637963(2) [2] 0.5924(1) [3] 

U A 0.5093833(3) [1] 0.5637961(3) [1] 1.2431(1)[1] 

U 6 0.5093834(3) [2] 0.5637962(3) [2] 1.7509(2) [7] 



Systematic errors are estimated by comparison with fits with 24 < L < 48, i.e., by evaluating 
(difference of the two fits)/(2 :r — 1) with x — 1/u + u « 2.3 for (3 C and x = uo 0.8 for 
R* (here we pessimistically assume that leading corrections dominate). The estimates of 
(3 C obtained by using different quantities are all consistent among each other. As our final 
result we take the one obtained from the data of Rz'- (3 C = 0.5093835(2) [3] for the 4 model 
at A = 2.07 and f3 c = 0.5637963(2) [2] for the ddXY model at D = 1.02, where the number 
in parentheses is the statistical error, while the number in brackets is the systematic error 
due to scaling corrections. In Table I we also report the estimates of the fixed-point values 
R* of the phenomenological couplings, which improve the results of Ref. 5. For the other 
values of A and D considered, we estimate f3 c by requiring that R(f3 C) L = 128) = R*; for R* 
we use the estimate of R* reported in Table I. Again, the best estimate is obtained from 
Rz; for this quantity scaling corrections are quite small. The results are reported in Table 
II, where the number in parentheses is the statistical error, while the number in brackets is 
the error due to the uncertainty on R*. 

D. FSS at fixed phenomenological coupling R 

Instead of performing the FSS analysis at fixed Hamiltonian parameters, one may analyze 
the data at a fixed value of a given phenomenological coupling R, as discussed in Sec. III. 
For this purpose we need to compute R(f3) in a neighborhood of f3 c . This could be done 
by reweighting the MC data obtained in a simulation at (3 = /3 run ~ (3 C . However, due to 
our enormous statistics, we could not store all results needed to reweight the data. Instead, 
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TABLE II: Estimates of j3 c for several values of A and D, from the FSS analysis of the MC 
data, and from the fit of MC data of x an d £ i n t ne HT phase using blAl approximants of the 
22nd-order HT series (MC-HT) of x and £, see Sec. VB. The results for A = 2.00 (0 4 model) and 
D = 0.90, 1.03 (ddXY model) are taken from Ref. 5. 



model 


FSS 


MC-HT from X 


MC-HT from £ 


</> 4 , A = 1.90 


0.5105799(4) [3] 






<p A , A = 2.00 


0.5099049(15) 






<f> 4 , A = 2.07 


0.5093835(2) [3] 






4 , A = 2.10 


0.5091503(3) [3] 


0.5091504(4) 


0.5091504(4) 


</> 4 , A = 2.20 


0.5083355(3) [4] 


0.5083361(4) 


0.5083363(4) 


ddXY, D = 0.90 


0.5764582(15)[9] 






ddXY,L> = 1.02 


0.5637963(2) [2] 


0.5637956(6) 


0.5637970(7) 


ddXY,L> = 1.03 


0.5627975(7) [7] 






ddXY,L> = 1.20 


0.5470376(17) [6] 


0.5470383(6) 


0.5470392(7) 


standard XY 


0.4541652(5) [6] a 







°This estimate is obtained from a fit with R Z (/3 C ,L) = R* z + cL~ u , with L min = 32, fixing R* z = 0.3203 
and uj = 0.785. This result is consistent with /3 C = 0.4541659(10) given in Ref. 35. 

we computed the first derivatives of R(/3) with respect to (3 and determined R(/3) by using 
its third-order Taylor expansion around /9 run . We checked that this is by far enough for our 
purpose. If we do not use the reweighting technique, it is enough to store bin-averages of 
the different quantities, significantly reducing the amount of needed disk space. Given R((3), 
one determines the value (3f such that R((3 — flf) — Rf. All interesting observables are then 
measured at their errors at fixed R = Rf are determined by a standard jackknife analysis. 
For compatibility with our previous study, 5 we choose Rzj = 0.3202 and R^j = 0.5925. 

This method has the advantage that it does not require a precise knowledge of the critical 
value (3 C . But there is another nice side effect: for some observables the statistical errors at 
fixed Rf are smaller than those at fixed (3 (close to f3 c ). This is due to cross-correlations and 
to a reduction of the effective autocorrelation times. For example, we find 

6rr[xk] -3.2, r 6r ^ ] ^4.5, (28) 



err[x|ij z=0 .3202] ""' err[x|i? e=0 .5925 
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err[C/ 4 |/3 c ] , n err^l^J 

: — ~ PS 1.9, : — ~ ~ 1.6 

err[f/ 4 | Rz=0 .3202j err[[/ 4 |.R 5 =0.5925] 

for the ddXY model, with a very small L dependence (within the last figure of the above- 
reported numbers). In the case of the 4 model we find slightly smaller improvements for 
the same quantities. We also mention that the gain is marginal for the derivatives of R 
considered in this paper. A reduction of the statistical errors when some quantities are 
measured at a fixed R has also been observed in other models. 36 



E. The leading correction-to-scaling exponent u 

In order to study corrections to scaling we analyze the value of a phenomenological 
coupling Ri at a fixed value of a second coupling R 2 . If /?/ is determined by R 2 (f3f) = R2J, 
we consider 

Ri = Ri((3 f ) (29) 

(the dependence on L is understood hereafter). Note that the large- L limit of R\ is uni- 
versal but depends on R 2 j- It differs from the critical value R\, unless R 2 j = R 2 - The 
phenomenological couplings that we consider are U4, Uq, Rz = Z a /Z pi R^ = £/L, and 
Ry = TL. 

In the 4> A model we define 

A(Ai,A 2 ) = ^(A 2 )- J R(Ai), (30) 

and analogously for the ddXY model, replacing A with D. Since 

R(\) = R* + C (A)L- W + . . . , 

we perform fits with the most simple ansatz 

A(Ai,A 2 ) = AcL- w , Ac = c(A 2 )-c(Ai), (32) 

and analogously for the ddXY model. A selection of such fits is given in Table III, where 
we report only results for C7 4 and R r at either a fixed value of Rz or R^. Using Uq instead 
of U4 leads to very similar results. R^ at a fixed value of R z is little useful, because the 
corrections to scaling and the statistical errors are relatively large. In order to get an idea of 
the corrections to scaling, we give results for the fit intervals 5 < L < 12 and 10 < L < 24. 
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(31) 



TABLE III: Fits to ansatz (32), using the data at Ai = 1.9 and A2 = 2.3 in the case of the (j) 4 
model and D\ = 0.9 and D2 = 1.2 in the case of the ddXY model. L m m and L max are the minimal 
and maximal lattice size that has been included in the fit. x 2 /d.o.f. is the chi square per degree of 
freedom of the fit. 



model 


Ri 


R2 


Lmin 




Ac 


to 


X 2 /d.o.f. 




u 4 


Rz 


5 


12 


-0.0209(2) 


0.825(4) 


0.8 








10 


24 


-0.0200(5) 


0.804(10) 


1.1 




u 4 


Rt 


5 


12 


-0.0210(2) 


0.775(4) 


0.5 








10 


24 


-0.0215(5) 


0.785(10) 


1.4 




Rr 


Rz 


5 


1 

12 


—0.0053(1) 


0.722(9) 


1.4 








10 


24 


-0.0060(5) 


0.775(37) 


0.7 


ddXY 


u 4 


Rz 


5 


12 


-0.0355(2) 


0.782(3) 


0.7 








10 


24 


-0.0349(6) 


0.775(7) 


0.9 




u 4 


Rt 


5 


12 


-0.0365(2) 


0.739(3) 


4.3 








10 


24 


-0.0386(7) 


0.764(7) 


0.6 




R r 


Rz 


5 


12 


-0.0099(1) 


0.708(5) 


3.6 








10 


24 


-0.0115(7) 


0.773(25) 


0.9 



The difference of the two results should be a rough estimate of the error due to next-to- 
leading corrections to scaling. As our final result we quote 

u = 0.785(20), (33) 

which includes (almost) all results for the interval 10 < L < 24. 

F. Determination of A* and D* 

To begin with, in Fig. 2 we show results for U 4 and R T at fixed R z = 0.3202, for various 
values of A in the 4 model, of D in the ddXY model, and in the standard XY model, vs. 

with u = 0.785. They show a clear evidence of the leading scaling corrections, and of 
the existence of optimal values A*, D* of A and D for which they are suppressed. We also 
note that -Rr is subject to larger next-to-leading scaling corrections; we shall return to this 
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FIG. 2: (above) and Rr (below) at fixed Rz = 0.3202 for various values of A (0 4 model), D 
(ddXY model) and the standard XY model, vs. with u = 0.785. 

point later. 

In order to determine A* and D*, we mainly use our data generated for A = 2.07 and 
D = 1.02. We fit them using various ansatze. The most simple one is 

R = R* + cL~ w , (34) 

where R is defined in Eq. (29). Eq. (34) includes only leading corrections to scaling; we fix 
uj = 0.785 as previously obtained. We may also include subleading corrections to scaling 
and consider 

R = R* + cL~ w + eL~ W2 . (35) 
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These fits are a bit problematic, since there are several scaling corrections that have similar 
exponents with 1.8<u n < 2.0, see Sec. III. In our fits we use both uo 2 = 1.8 or uo 2 = 2. 

We first perform fits of types (34) and (35) for the two models 4 and ddXY separately. 
As R we consider at fixed Rz, at fixed R^, Uq at fixed Rz, U% at fixed R^, and _R-r at 
fixed Rz- The results for R* are, as required by universality, consistent for the two models. 
Hence we take our final results from joint fits of the results for both models. For instance, 
in a joint fit with ansatz (34) there are three free parameters: R* and a correction-to-scaling 
amplitude for each of the 4 and ddXY models. 

In order to determine A* (and analogously D*) we assume c(A) to be linear in a neigh- 
borhood of A* and write 

c(A) « Ci(A- A*), (36) 

so that 

A* = A - -c(A). (37) 

Ci 

We use A = 2.07, the value for which we have most of the simulations, and determine Ci by 
using 

dc 



Cl ~d\ 



c(A = 2.3) - c(A = L9) _ 



2.3-1.9 



In the ddXY model we use the same formulas with D — 1.02 and 



dc 

Cl ~dD 



c{D= 1.2) -c{D=^ 



1.2 0.9 



In order to determine the needed values of c(A) and c(-D), we fix uo = 0.785 (and, to estimate 
errors, uo = 0.765, uo = 0.805). The results of the fits with ansatz (32) for uo = 0.785 are 
summarized in Table IV. The final estimate is taken from the fits with 12 < L < 24. The 
comparison with the fits with 6 < L < 12 gives us an idea of the error due to subleading 
corrections. It is small enough to be ignored in the following. 

We also checked whether the linear approximation is sufficiently accurate by determining 
the derivatives (38) and (39) from other pairs of values of A and D. The error due to the 
linear extrapolation is approximately 10%, which is negligible for the purpose of determining 
A* and D*. In Fig. 3 we plot the results for A* and D* as functions of the L m i n used in the 
fits to (34) and (35). The results of the fits to Eq. (34) show a systematic drift and become 
stable only for L min = 30. This systematic variation is mostly due to the next-to-leading 
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TABLE IV: Fits to ansatz (32) 


as in Table III, but keeping oj fixed. 


Ac = c(A = 2.3) 


„(\ 1 Q\ 


for the model and Ac = c(D 


= 1.2) - 


c(D = 0.9) for the ddXY model. 




rnnnpl ??i 

lllUVlvl J. L I 


R2 


UJ 


L ■ 

-^min 


L 


Ac 


X 2 /cLo.f. 


6 4 TL 




n 78^ 

u. 1 00 


6 


12 


—0 01925(31 


5.8 








1 9 


94 


—0 01 898(51 


1.5 






u.ouo 


p. 

yj 


12 


—0 02006(31 


1.5 










94 


— 09009(51 


1.2 






n 78^ 




12 


—0 01847(31 


13.6 








1 9 

.1. .1. 


94 


—0 01 799(51 


2.4 


U4 


T?? 


U. 1 00 




1 9 

_L Z. 


—0 091 45(31 

\J .\JZj ± L ±tJ IOI 


0.2 








12 


24 


—0 091 45(51 


1.6 


TJr 


Ry 


n 78^ 


G 


12 


—0 0688(11 


7.7 








1 9 


94 


—0 0fi78(9l 


1.5 


Ur 




n 78^ 


p. 

yj 


12 


—0 0763(1 1 


0.2 








1 9 


94 


— D D76^(9l 


1.6 






n 78^ 

U. 1 00 


a 


12 


—0 00608(91 


2.5 








1 9 


9/1 


— U.UUDIO^O ) 


0.8 


ddXY TL 


HZ 


U. (80 


u 


12 


—0 03576(31 


0.8 








1 9 
1Z 


94 


— W.UOO ( oyoj 


0.9 






U.oUD 




12 


— 017270) 


6.5 








1 9 


94 


— U.UO / l ^0 J 


3.4 






U. ( DO 




12 


U.UL>t:<J1 1 U J 


4.4 








1 9 

J. ^ 


94 


—0 03393(71 


0.2 


TL 


u < 


n 78^ 




12 


—0 04021 (41 


12.7 








1 9 


94 


—0 04077(91 


1.0 


Ur. 


Rz 


n 78^ 

U. 1 00 


6 


12 


-0.1271(1) 


1.7 








12 


24 


-0.1276(3) 


1.2 


Ur 


Rt 


0.785 


6 


12 


-0.1423(1) 


22.3 








12 


24 


-0.1447(3) 


1.3 


Ry 


Rz 


0.785 


6 


12 


-0.01157(2) 


7.2 








12 


24 


-0.01177(8) 


1.2 
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FIG. 3: Determination of A* and D*: (a) results from fits of U4 at Rz = 0.3202 to a + cL u with 
uj = 0.785; (b) results from fits of U4 at R^ = 0.5925 to a + cL~ w with lu = 0.785; (c) results from 
fits of U 4 at R z = 0.3202 to a + cL~ u + eL" W2 with u = 0.785, uj 2 = 1.8; (d) results from fits of U A 
at = 0.5925 to a + eL -£J + eL~ U2 with u> = 0.785, U2 = 1-8. The dashed lines indicate our final 
estimates. 
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corrections and indeed, fits to Eq. (35) are less dependent on L min and give fully consistent 
results. As our final result we quote 

A* = 2.15(5), D* = 1.06(2), (40) 

which correspond to 

(^1^=0.3202)* = 1.24281(10), (C7 4 lfl^o.5925)* = 1.24277(10). (41) 

Consistent results for A* and D* are obtained by analyzing Uq at Rz = 0.3202 and = 
0.5925. Note that the estimates of A* and D* are slightly larger than those obtained in our 
previous work, 5 where we reported A* = 2.07(5) and D* = 1.02(2). Since the larger statistics 
and larger lattice sizes allow us to achieve a better control of all sources of systematic errors, 
and in particular of the next-to-leading scaling corrections, we are confident that the new 
estimates (40) are now correct with the quoted errors (that are as large as those reported in 
our previous work). 

We also performed some MC simulations of the standard XY model up to lattice size 
L = 96. Using the estimates for uj and Ul obtained above, one-parameter fits of the MC 
data give (taking data for 24 < L < 96). 

^1^=0.3202 = 1-24281 - 0.1014(4)L-°- 785 (42) 

and 

£4^=0.5925 = 1-24277 - 0.1138(4)L-°- 785 . (43) 

We can use these results to obtain a conservative upper bound on the ratios |c(A = 
2.15)/c(XY)| and \c(D = 1.06)/c(XY)| that are independent on the quantity one is con- 
sidering. Using the estimate of <9c/<9A|a=a* and taking into account that the error on A* 
is ±0.05, we infer that the leading scaling-correction amplitude of U4 at Rz = 0.3202 for 
A = 2.15 satisfies |c(A = 2.15) | < 0.0024. The same bound is obtained in the ddXY model 
for c(D = 1.06). This implies |c(A = 2.15)/c(XY)| < 0.0024/0.1014 w 1/42 and an anal- 
ogous bound for \c(D = 1.06)/c(XY)|. A similar calculation also shows that |c(A = 2.07) | 
and \c(D = 1.02) | are at least 20 times smaller than |c(XY)|. 
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G. Next-to-leading scaling corrections 



In this subsection we present evidence of next-to-leading scaling corrections characterized 
by an exponent cj mo ~ 2, as expected due to the presence of several irrelevant perturbations 
with y ~ —2, as discussed in Sec. III. In particular, this provides a robust evidence of the 
absence of 1/L analytic corrections, see discussion at the end of Sec. III. 

We first construct improved variables Ua\r z =0.3202, £41^=0.5925, and -Rt|,r z =o.3202, which 
do not have leading scaling corrections. In the ddXY model (analogous formulas hold in the 
4 model by replacing D with A) we consider 

Amp ~ R(D 1 ) X R(D 2 ) 1 - X . (44) 

Expanding R as in Eq. (34) and using a linear approximation for c(D), c(D) ps c±(D — D*), 
as in the previous section, we obtain for L — > 00 



-Rimp ~ R 



l + ^(xD 1 + (l-x)D 2 -D*) + ... . (15) 



Thus, if we take x = x* = (D* — D 2 )/(D 1 — D 2 ), the leading scaling correction cancels. In 
the ddXY model we use the data at D\ = 1.02 and D 2 = 1.2, while in the 4 model we 
combine our data for A x = 1.9 and A 2 = 2.3. 

In Fig. 4 we show the improved combination (44) in the case of the Binder cumulant U4 
as a function of L~ 2 . For t7 4 | ^=0.5925 we observe straight lines to a good approximation. On 
the other hand, for £/4|i? z =o.3202 there is a clear bending of the curves, indicating that, in the 
given range of lattice sizes, corrections with exponent u/ > 2 are contributing significantly. 
We should note that these results do not provide a completely independent information, 
since we have already used oo 2 ~ 2 as input in the determination of D* and A* in the 
previous section. Therefore, they only provide a consistency check. This is not the case for 
R r = TL, since it was not used for the determination of D* and A*. In Fig. 4, -Rt|r z =o.3202 
is plotted as a function of L~ 2 . We clearly observe a straight line. Hence corrections with 
oj 2 pa 2 clearly dominate in the whole range of lattice sizes that are shown. There are 
no leading corrections to scaling and, as expected from general RG arguments, also no 
corrections proportional to L^ 1 . If there were a correction proportional to the ratio of 
its amplitude with the amplitude of the correction proportional to (i.e., the leading one) 
would not be the same in different quantities. In other words, if corrections proportional to 
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FIG. 4: Estimates of the improved quantity (44) in the case of the Binder cumulant U4 at 
R z = 0.3202 and = 0.5925 (above), and of R r at R z = 0.3202 (below), vs. L~ 2 . 

L^ 1 and effectively cancel for U4 at our estimates of D*, A* and for the range of values 
of L considered, there is no reason why this should also happen for R-j. Hence we conclude 
that our numerical results confirm the theoretical argument that no L~ l corrections are 
present in FSS for periodic boundary conditions. 
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TABLE V: Fits to ansatz (49) of Sjj 4 at Rz = 0.3202 for the ddXY model, where we have used 
D\ = 0.9 and D 2 = 1-2. We fix to = 0.785. L m \ n and L max are the minimal and maximal lattice 
size included in the fit. 



min 


-^max 


a(D = 1.2)/a(D = 0.9) 


c(D = 1.2) - c{D = 0.9) 


xVcLo.f. 


6 


12 


0.9811(4) 


-0.0735(18) 


1.4 


6 


24 


0.9809(2) 


-0.0723(12) 


1.0 


8 


24 


0.9806(3) 


-0.0699(21) 


0.9 


10 


24 


0.9803(5) 


-0.0677(33) 


1.1 


12 


24 


0.9803(7) 


-0.0672(54) 


0.3 



H. The critical exponent v 

Here, we compute the critical exponent v from a FSS analysis of the derivative S\ = 
dRi/df3 at a fixed value of another (or the same) phenomenological coupling R 2 . 

For this purpose we define an improved quantity that does not have leading scaling 
corrections. Since S 1 behaves as 

Si = a{\)L l ' u [l + c(A)L- w + . . .] (46) 

for L — > oo, if we take A A*, so that c(A) is small and the linear approximation c(A) ~ 
Ci(A — A*) works well, an improved variable is simply 

Si )imp (A) » 51(A) [1 - Ci(A - A*)L-] . (47) 

We compute Ci using 



« C(A1) ~ C(A2) , (48) 

A=A* A l ~ A 2 

where Ai and A2 are sufficiently close to A* so that the linear approximation works well. We 



dc 

Cl= dX 



estimate the difference c(Ai) — c(A 2 ) from fits of 

Agi(Ax, A 2 ) ^ S f ^^ = a ~^\{l + [c(A 2 ) - c(A0] L~ w + ...}. (49) 

In the 4 model we take Ai = 1.9, A 2 = 2.3, and A = 2.07 and fix u — 0.785. The same 
formulas hold in the ddXY model: we take D 1 = 0.9, D 2 = 1.2, D = 1.02. 

As an example, let us discuss the determination of c(D = 0.9) — c(D = 1.2) for S*^ 
at i?z = 0.3202. In Table V we report some results of fits with ansatz (49). They are 
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quite stable when L min and L max are varied. This indicates that leading scaling corrections 
dominate in the difference and that subleading corrections vary little with D. As our final 
result we take 

c(D = 1.2) - c(D = 0.9) = -0.067(7), (50) 

which should take into account both statistical and systematic errors. The error due to 
the uncertainty on u as well as the error due to the linear approximation are negligible. 
We repeat this procedure for all quantities of interest. Typically, the amplitude differences 
analogous to (50) can be determined with an error of approximately 10% for the ddXY 
model and 15% for the 4 model. 

To determine v, we fit the data of S\ = Su 4 at R 2 = Rz = 0.3202 to 

dR 1 



Si 
and 



R2 =const 



d(3 



= aL 1 ^, (51) 

i?2 =const 



^il« 2 =co„st = ^(l + ^ 2 ) (52) 

where we consider either 0J2 = 1.8 or u} 2 = 2. Assuming that there are no leading scaling 
corrections, the fits of the original data at A = 2.07 and D = 1.02 give the result v = 
0.67181(12) and v = 0.67195(13) for A = 2.07 and D = 1.02, respectively. The errors take 
into account the results obtained by using the two ansatze and the L min dependence. 

In order to evaluate the effect of the residual leading corrections to scaling, we repeat the 
same fits for S^mp computed by using Eq. (47). We obtain estimates of v that are smaller 
by roughly 0.0002. This change depends slightly on the ansatz and L mm . These results 
allow us to give an effective estimate of v as a function of D and A, the A and D dependence 
being due to the residual leading corrections to scaling that are not taken into account in 
the fit. Since the estimates of the improved quantities correspond approximately to those 
that would be obtained by using data at A = A* or D = D*, we obtain 

v = 0.67181(12)- 0.0022(3) x (A -2.07) for the 4 model, (53) 
v = 0.67195(13) - 0.0061(9) x (D - 1.02) for the ddXY model. (54) 

The A and D dependence is that corresponding to the fit without corrections with L min = 
16. The error on the linear coefficient is due to the error on dc/d\. Using the estimates 
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A* = 2.15(5) and D* = 1.06(2), we obtain the results 



v = 0.67163(12) [11] for the 4 model, 
v = 0.67171(13) [12] for the ddXY model, 



(55) 
(56) 



where the first error is statistical and the second one is due to the uncertainty on A* and 
D*. The analysis of S Ua at = 0.5925 gives analogous results. Other quantities provide 
consistent, but less precise, estimates. 

Finally, we analyze the derivative Sr t = dRr/d(3 of the helicity modulus. In this case 
subleading corrections are smaller when considering the data at (3 — /3 C , rather than at a 
fixed phenomenological coupling. Since we have little data for the 4 model at A = 2.07, 
we only give results for the ddXY model. We first compute the improved slopes Sij Ti i mp 
by using Eq. (47). Fits to ansatz S = aL 1 /" have a small xVd.o.f. (d.o.f. is the number 
of degrees of freedom of the fit) already from rather small L min , indicating that subleading 
corrections are rather small in this quantity. Fitting the data with 16 < L < 128, we 
obtain v = 0.67200(15) at D = 1.02 and v = 0.67173(15) at D = 1.06 (more precisely, for 
the improved slope). Therefore, taking also into account the error on D*, we might quote 
v = 0.6717(3) as final result of this analysis, which agrees with the results (55) and (56). 
We have also checked the dependence of this result on the estimate of f3 c , finding that the 
error due the uncertainty of f3 c is definitely smaller than the error on v quoted above. 

I. Eliminating leading scaling corrections from the derivative of phenomenological 
couplings 

We now consider a combination of the derivative Si = dRi/dfi of two phenomenological 
couplings Ri, at (5 C or at fixed Rj, 



and show that one can choose a value p such that this quantity is improved — no leading 
corrections to scaling — for any A orD. Moreover, the computation of p does not rely on any 
estimate of u. 

For L — > oo, Si at (3 C or at a fixed Rj, behaves as 



Si{\)\ 



S 2 (A)| 



i-P 



(57) 



Si(A) 



CLi(\)L 



(58) 
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Therefore, we have 



|^(A)r \S 2 (\)\ 1 ~ p = \ ai (X)\ p |a 2 (A)| 1_p L}l u {1 + [p Cl (A) + (1 -p)c 2 (A)] L~» + . . .} .(59) 



An improved quantity is obtained by taking p = p*, where 



p* Cl (A) + (l-p*)c 2 (A) = 0. 



(60) 



Note that, since the ratios ci(A)/c2(A) are universal, i.e., independent of A (or D), the 
optimal value p* is universal. 

Now we show howp* can be accurately computed. We consider ratios of Si (A) at different 
values of A: 

Si(X 2 ) aj(A 2 ) 



{l + [ Ci (A 2 )- Ci (A 1 )]L- w + ...} 
Si{\i) aj(Ai) 

Due to the universality of the amplitude ratio we have 

c 2 (Ai) c 2 (A 2 ) c 2 (A 2 ) - c 2 (Ai) 



ci(Ai) ci(A 2 ) ci(A 2 ) - ci(Ai) ' 



(61) 



(62) 



Therefore, 



Si(A 2 



S 2 (X 2 ) 
S 2 (X 1 ) 



i-p 



«i(A 2 ) 
ai(Ai) 



Qi(A 2 ) 
ai(Ai) 



i-p 



(63) 



x {1 + [p ( Cl (A 2 ) - ci(Ai)) + (1 - p) (c 2 (A 2 ) - c 2 (Ax))] L" w + . . .} 



We can obtain the desired value of p* by imposing that the combination 



Si(X 2 ) 
Si(A0 



S 2 (X 2 ) 
S 2 (A0 



i-p 



(64) 



is L independent. This procedure does not require any knowledge of uj and only assumes 
that leading scaling corrections dominate in the considered range of lattice sizes. 
The optimal pair of slopes Si turns out to be 



Si 



dp 



So 



8Ua 



0(3 



(65) 



i? z =0.3202 



R z =0.3202 

essentially because the amplitudes of their leading scaling correction have opposite signs. In 
order to determine the corresponding value of p*, we perform fits with ansatz (64) using the 
data at A x = 1.9 and A 2 = 2.3 for the 4 model, and at D x = 0.9 and D 2 = 1.2 for the ddXY 
model. Some results are reported in Table VI. The estimates of p* from the two models are 
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FIG. 5: Results for the exponent v obtained by fits for several values of L m j n : (a) to aL x l v and 
(b) to aL l l v {\ + eL-^) with uo 2 = 1-8. 2? 



TABLE VI: Fits to ansatz (64) of the improved combination of Sjj 4 and Sr z at Rz = 0.3202. L m [ n 
and L max are the minimal and maximal lattice size included in the fit. For a detailed discussion, 



see the text. 

model L min L max p X 2 /d.o.f. 

4 6 12 0.703(11) 0.93 

8 24 0.718(11) 1.37 

12 24 0.744(30) 1.61 

ddXY 6 12 0.736(6) 1.87 

8 24 0.720(8) 1.29 

12 24 0.702(19) 0.02 



consistent, as required by universality. We take p* = 0.72(3) as our final estimate (from fits 
with L min = 8, while the error is estimated by varying the fit range). 

Using this estimate for p*, we construct the improved combinations (59) at A = 2.07 for 
the 4 model and at D = 1.02 for the ddXY model. Since these values of A and D are close 
to the optimal values A* and D*, leading scaling corrections are quite small. Therefore, 
the uncertainty on p* is negligible with respect to the final error of our estimate for v. In 
Fig. 5 we show results for the exponent v obtained by fits to the functions aL 1 ^ and to 
aV~l v {\ + eL^' 2 ) with u 2 = 1.8, for several values of L min . Guided by Fig. 5, we take 

v = 0.6718(2) for the 4 model, (66) 
v = 0.6717(3) for the ddXY model (67) 

as our final estimates. Moreover, combined fits applied to both 4 and ddXY models give 
the estimate v = 0.6718(2). 

As further check, we apply this method to the standard XY model, where leading cor- 
rections to scaling are large. In Fig. 6 we show results for the critical exponent v from MC 
simulations of the standard XY model up to L = 96, as obtained by fits to the simple ansatz 
S = aL x l v of the data of Sr z at Rz = 0.3202, Su 4 at Rz = 0.3202, and their combination 
(59) at p = p* = 0.72. The first two sets of results clearly disagree with the estimate of v 
from improved Hamiltonians and also between themselves (thus, they are inconsistent). In- 
stead, the analysis of their improved combination provides perfectly consistent results, giving 
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FIG. 6: Results for the critical exponent v from MC simulations of the standard XY model, as 
obtained by fits ot the data of the derivative of Rz at Rz = 0.3202 (a), the derivative of C/4 at 
Rz = 0.3202 (b), and by using the method that combines them to eliminate the residual leading 
scaling corrections (c). The two dashed lines correspond to our final estimate v = 0.6717(1). 

further support for the validity of the method and confirming the accurate determination of 
p*. 
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V 
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J. The critical exponent v from the finite-size scaling of the energy density 

We also derive estimates of v from the FSS of the energy density. On finite lattices of 
size L 3 , the free energy density behaves as 

/(A L) = f ns (f3) + f s ((3, L), (68) 

where the nonsingular part of the free energy density f ns does not depend on L, apart from 
exponentially small contributions. The singular part is expected to behave as 

f s ((3,L) = L~ d g(tL 1 /»,u 3 L-",...). (69) 

The energy density is obtained by taking the derivative with respect to (3, 

E(/3, L) = E ns {(3) + L- d+1 l v g e (tL l '\ u 3 L^, ...) + ... (70) 
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Setting (3 = (3 C , we obtain for L — > oo the expansion 

E(/3 C , L) = E ns {f3 c ) + aL~ d+1 ' v {l + cL" w ) + . . . (71) 

In Fig. 7 we show the results of the fits to Eq. (71) without correction terms (i.e., with c = 0) 
for the 4 and ddXY models respectively at A = 2.07 and D = 1.02. Our final estimates of 
v from the scaling of the energy density (obtained from fits of data for 12 < L < 128) are 

v = 0.6717(2) for the 4 model, (72) 
v = 0.6715(3) for the ddXY model. (73) 

Our data at A = 1.9,2.3 in the case of the 4 model and D = 0.9,1.2 for the ddXY 
model are not sufficient to get a reliable estimate of the systematic error due to the leading 
scaling corrections. However, we can estimate it by a comparison with the standard XY 
model. Fitting the data of the standard XY model to the simplest ansatz without scaling 
corrections (using the value of f3 c reported in Table II and data for 8 < L < 96), we obtain 
v = 0.6701(2), i.e., the exponent is underestimated by approximately 0.0015. We may use 
this difference to estimate the leading scaling corrections amplitude in the standard XY 
model. Taking into account that the largest lattice for the XY model has L = 96 instead 
of L = 128, we obtain c « 0.0015 x (128/96)" - 785 = 0.0012, where we have also taken into 
account the difference of the sizes L used in the fits. Since the leading correction amplitudes 
at D — 1.02 and A = 2.07 should be smaller by a factor of approximately 20 than in the 
standard XY model, we may have a shift by 0.00006 in our estimates (72) and (73) for v 
due to leading corrections, which is much smaller than their errors. 

Note that, in the case of the energy, analyses at fixed phenomeno logical coupling do not 
work. Since 

R(L, (3)=R* + a(P - p c )L l ' v + bL' w + . . . , (74) 

and therefore 

R f = R* + a (p f _ [3 c )L 1/u + bL-" + . . . , (75) 

then 

[3 f = [3 c + a-\R f - R*)L- 1/u - (b/a)L-^ v - u + ... (76) 

Plugging Pf in the nonsingular part of the energy, one obtains a term 0(L _1 / 1/ ), which, in 
this particular case, completely masks the interesting scaling term, since \jv — 3 ~ —1.51 
and ps -1.49. 
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FIG. 7: Results of fits for the FSS of the energy for the <j> model at A = 2.07 and the ddXY model 
at D = 1.02. 

K. The critical exponent ij 

We compute the critical exponent r? from the FSS behavior of the magnetic susceptibility 
X either at Rz = 0.3202 or at = 0.5925. As in the analyses of the derivatives Sr, we first 
compute an improved quantity for x that does not have leading scaling corrections. Since 



X = aL 2 - r > (1 + cL-» + ...) 



(77) 



close to the improved value A* where c is small and a linear approximation c(A) = Ci(A — A*) 
suffices, we can take 



Ximp = x(A) (1 - c(X)L-) « X (A) [1 - Cl (A - \*)L-"] . 



(78) 



To compute Ci we consider the ratio 

X(^2) 1^=0.3202 0(A 2 ) 



{l + [c(A 2 )-c(Ai)]L- w + ...} 



X(^l) lfl z =0.3202 

with A 2 = 2.3 and A2 = 1.9. Fits to Eq. (79) allow us to estimate the difference c(A 
2.3) -c(A = 1.9). Then 



(79) 



dc 

Cl= d\ 



A=A* 



c(A = 2.3) -c(A = 1.9) 
2.3 - 1.9 ' 



(80) 



Analogous equations can be written for the ddXY model and at fixed R^. 
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FIG. 8: Results for the exponent r\ obtained by fits of Ximp f° r several values of L mm : (a) to 



aL 1 ^ and (b) to aL x l v {\ + eL -Wa ) with u 2 
7} = 0.0381(3). 



1.8. The dashed lines indicate our final FSS estimate 



Then, we fit the data for the improved x [ we use ^ = 2.07 in Eq. (78)] using the ansatze 



X 



aL 2 -\ 



lm Pl_R z =0.3202 

Ximp| fiz=03 202 = a L v (l + eL 2 ) , 



(81) 
(82) 



where we fix either 0J2 = 1.8 or 002 = 2. In Fig. 8 we show the results of the fits of 
Xim P |R z =o.3202 for the 4 model vs. the minimum value L min of L allowed in the fits. Similar 
results are obtained for the ddXY model, using D = 1.02 in Eq. (78). In all cases, the fits 
with ansatz (81) have a large x 2 /d-o.f. for L min < 40. Moreover, the resulting values of rj 
appear to slightly increase with increasing L min . In contrast, the fits allowing for subleading 
corrections are more stable and give a x 2 /d.o.f. close to 1 already for L min > 10. We also 
mention that fixing instead of Rz gives slightly lower values for rj, by about 0.00005 for 
the fit to (81), and 0.00015 for the fits to (82). A comparison of A = 2.07 and A = 2.15 
indicates that a possible error due to the uncertainty of A* should be approximately 0.00005. 
Our final FSS estimate obtained from both <b A and ddXY models is 



77 = 0.0381(3). 



(83) 
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V. CRITICAL EXPONENTS FROM IMPROVED HIGH-TEMPERATURE EX- 
PANSIONS 



A. High-temperature expansions 

High-temperature (HT) expansion in powers of the inverse temperature (3 is one of the 
most effective lattice approaches to investigate the critical behavior in the HT phase. We 
consider a general class of models defined on a simple cubic lattice by the Hamiltonian 

(xy) x 

where f3 = 1/T, (xy) indicates nearest-neighbor sites, <p x = (<f)x ,<f>x ) is a two-component 
real variable, and V{<p 2 ) is a generic potential satisfying appropriate stability constraints. 

Using the linked-cluster expansion technique (see Refs. 5,37 for details), we extended the 
HT computations of Ref. 5 by adding a few terms in the HT series. The two main steps 
of the algorithm are the generation of the graphs and the evaluation of the contribution of 
each graph. The first step is limited by memory, and was run on a computer with 16 Gbytes 
of RAM. The second step is limited by processing time; it was parallelized and required 
approximately 4 years of CPU-time on a single 2.0 GHz Opteron processor. We computed 
the 22nd-order HT expansion of the magnetic susceptibility and the second moment of the 
two-point function 

X = ]T(0 a (O)0 a (x)>, (85) 

X 

m 2 = ^2x 2 ((f) a (0)(j) a (x)), 

X 

and therefore, of the second-moment correlation length 

e = p. (86) 

Moreover, we computed the HT expansion of some zero-momentum connected 2j-point 
Green's functions X2j 

X2j= (^i(°)^«i( a; 2)--0a J (a;2i-i)0a J (a;2j)) c (87) 

X2, — ,X2j 

(x — Xz)- We computed %4 to 20th order, xe and %8 to 18th order. 

The HT series for the general model (84) will be reported elsewhere. In the following we 
will restrict ourselves to the 4 and ddXY models, cf. Eqs. (9) and (10). 
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B. Critical exponents from the analysis of the HT expansion of improved models 

In our analysis of the HT series we consider quasi-diagonal first-, second-, and third- 
order integral approximants (IAl's, IA2's and IA3's respectively), and in particular biased 
IAn's (blAn's) using the most precise available estimate of (3 C . The FSS estimates of f3 c are 
reported in Table II. We refer to Ref. 5, and in particular to its Appendix B, for details on 
the HT analysis and the precise definition of the various integral approximants. A review 
of methods for the analysis of HT series can be found in Ref. 38. 

The leading nonanalytic corrections are the dominant source of systematic errors in HT 
studies. Indeed, nonanalytic corrections introduce large and dangerously undetectable sys- 
tematic deviations in the results of the analysis. Integral approximants can in principle cope 
with an asymptotic behavior of the form (13); however, in practice, they are not very effec- 
tive when applied to the series of moderate length available today. As shown in Refs. 5,21,23, 
analyses of the HT series for the improved models lead to a significant improvement in the 
estimates of the critical exponents and of other infinite- volume HT quantities. The crux of 
the method is a precise determination of the improved value of the parameter appearing in 
the Hamiltonian. In this respect FSS techniques appear quite effective, as we have shown in 
the preceding section. A further improvement is achieved by biasing the HT analysis using 
the available estimates of f3 c . 

Our working hypothesis is that, with the series of current length, the systematic errors, 
i.e., the systematic deviations that are not taken into account in the HT analysis, are 
largely due to the leading nonanalytic corrections, especially when they are characterized 
by a relatively small exponent, as is the case in the 3-d XY universality class where A = 
ojv ps 0.53. Therefore, improved models are expected to give results with smaller and, 
more importantly, reliable error estimates. The systematic errors in our analyses are related 
either to next-to-leading nonanalytic scaling corrections or to our approximate knowledge 
of A* and D*. Since A 2 = uj 2 v ~ 1-2, we will assume that next-to-leading corrections 
do not play much role, and we will only take into account the residual leading corrections 
proportional to A — A* (D — D* for the ddXY model). Of course, this hypothesis requires 
stringent checks. A nontrivial check is achieved by comparing the results obtained using 
different improved models: If the hypothesis is correct, they should agree within error bars. 

As an additional check of our results we compare IAl's of the 22nd-order HT series of 
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X and £ 2 with high-statistics MC results. We simulated the <fi 4 model at A = 2.1, 2.2 and 
the dd-XY model at D — 1.02, 1.2 in the HT phase (f3 < (3 C ). We alternate single-cluster 
updates and local Metropolis and overrelaxation updates. In order to obtain negligible finite- 
size effects (i.e., orders of magnitude smaller than the statistical errors) we used lattices of 
size L > 10£ throughout. We obtained infinite- volume estimates up to £ ~ 30 on a 350 3 
lattice. 

In Fig. 9 we show MC data for the <p A model at A = 2.1 from (5 = 0.493, where f = 
4.1825(2), to the largest (3 value (3 = 0.5083, where f = 30.453(10). blAl's using the FSS 
estimate (3 C = 0.5091503(5) provide perfectly consistent results, for example £ = 30.449(1) [7], 
where the first error is related to the spread of the approximants and the second one to the 
uncertainty on f3 c . We also obtain an independent estimate of f3 c by fitting the MC data of 
£ to blAl's with f3 c taken as a free parameter. The resulting estimate (3 C = 0.5091504(4) 
(with xV d -°- f ~ °-9) is perfectly consistent with the FSS estimate (3 C = 0.5091503(5). The 
corresponding curve is drawn in Fig. 9. An identical result, i.e., (3 C = 0.5091504(4), is found 
by fitting the MC data of x, which shows that the agreement with the FSS analysis is not 
just by chance. We performed similar analyses also for A = 2.2, and in the case of the 
ddXY model for D = 1.02 and D = 1.20. The results denoted by MC+HT are reported 
in Table II. The comparison with the FSS estimates is overall satisfactory. This successful 
analysis should be contrasted with the case of the standard XY model, where the fit of the 
MC results in the HT phase by blAl does not provide acceptable results: indeed, most of the 
approximants are defective. This fact may be explained by the presence of sizable leading 
scaling corrections. 

In order to determine the critical exponents 7 and u, we analyze the 22nd-order series 
of x an d 21st-order series of / f3 respectively, using blAn's with n = 1,2,3 biased at the 
best available estimate of (5 C . Estimates of the exponent 77 can be obtained by using the 
scaling relation r\ = 2 — 7///. More precise estimates of the product r\v can be obtained by 
using the so-called critical-point renormalization method (CPRM) applied to the series of 
X and £ 2 / '(3, see Ref. 5 for details. For the 4 and ddXY models we performed analyses at 
different values of the A and D to determine the dependence of the effective exponents on 
A and D due to the residual leading corrections to scaling that are not taken into account 
in the analysis. In Table VII we report some intermediate results, i.e., the results for each 
blAn analysis, which will then lead us to our final estimates reported below. We closely 
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FIG. 9: High-temperature MC data for the correlation length of the </> 4 model at A = 2.1, and the 
result of the fit using blAl of the 22nd-order HT series of £ 2 . 

follow Ref. 5, so we refer to it for details. We only mention that the set of blAra's that we 
consider are those with q = 2, s = 1/2, n a = 2 in the definitions reported in Ref. 5. 

In the case of the 4 model, using the results of the blAn analysis at A = 2.10 and 
A = 2.20, and assuming a linear dependence on A in between, we obtain the IHT results 

v = 0.67161(4)[2] + 0.0018(A - 2.10), (88) 
7 = 1.31768(3)[5] + 0.0021(A - 2.10), 
r]v = 0.02556(5) + 0.0013(A - 2.10). 

The central value at A = 2.1 is taken from the bIA2 and bIA3 analyses, the number in 
parentheses is basically the spread of the approximants at A = 2.10 using the central value 
of p c (we use (3 C = 0.5091504(4)), while the number in brackets gives the systematic error due 
to the uncertainty on (3 C . The dependence of the results on the chosen value of A is estimated 
by assuming a linear dependence, and evaluating the coefficient from the results for A = 2.2 
and A = 2.1, i.e., from the ratio [Q(X = 2.2) — Q(X = 2.1)] /0.1, where Q represents the 
quantity at hand. Consistent results are obtained by using the pair of values 2.07, 2.20 or 
1.90, 2.10 instead of 2.10, 2.20. In the case of the ddXY model, using the results of the blAn 
analysis at D — 1.02 and D = 1.20, we obtain 

v = 0.67145(6) [2] + 0.0046(£> - 1.02), (89) 
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TABLE VII: Results for the exponents 7, v, and rj as obtained by using blAn. 



model 




approximants 


7 




r)v from CPRM 


<p A , A = 1.90 


0.5105799(7) 


blAl 


1.31718(3) [10] 


0.67116(10)[3] 


0.02531(4) 






bIA2 


1.31726(8) [12] 


0.67119(5)[3] 


0.02527(10) 






bIA3 


1.31723(3) [11] 


0.67118(14) [4] 


0.02528(5) 


4> 4 , A = 2.07 


0.5093835(5) 


blAl 


1.31757(2)[7] 


0.67153(4) [3] 


0.02553(5) 






bIA2 


1.31759(4)[7] 


0.67155(6)[3] 


0.02552(12) 






bIA3 


1.31758(2)[7] 


0.67153(5)[3] 


0.02552(4) 


c/) 4 , A = 2.10 


0.5091504(4) 


blAl 


1.31767(3)[6] 


0.67160(3)[2] 


0.02557(5) 






bIA2 


1.31769(5)[5] 


0.67163(6) [2] 


0.02556(12) 






bIA3 


1.31767(3)[5] 


0.67160(4) [2] 


0.02556(5) 


c/> 4 , A = 2.20 


0.5083355(7) 


blAl 


1.31787(2)[8] 


0.67178(2)[5] 


0.02569(6) 






bIA2 


1.31789(5) [10] 


0.67181(7)[5] 


0.02567(12) 






bIA3 


1.31788(5)[9] 


0.67179(4) [4] 


0.02569(5) 


ddXY,L> = 1.02 


0.5637963(4) 


blAl 


1.31757(15)[5] 


0.67141(4)[3] 


0.02547(5) 






bIA2 


1.31746(9) [4] 


0.67143(6) [2] 


0.02550(23) 






bIA3 


1.31745(12)[5] 


0.67149(9) [3] 


0.02533(19) 


ddXY,L> = 1.03 


0.5627975(13) 


blAl 


1.31767(15)[15] 


0.67148(4) [7] 


0.02550(4) 






bIA2 


1.31756(9) [14] 


0.67150(5) [6] 


0.02552(23) 






bIA3 


1.31755(10)[15] 


0.67155(8)[9] 


0.02537(18) 


ddXY,L> = 1.20 


0.5470388(11) 


blAl 


1.31871(5) [14] 


0.67227(8) [7] 


0.02602(4) 






bIA2 


1.31872(11)[13] 


0.67232(10) [6] 


0.02597(28) 






bIA3 


1.31868(13)[13] 


0.67224(10) [8] 


0.02575(52) 
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7 = 1.31746(11)[5] + 0.0070(L> - 1.02), 
•qv = 0.02547(14) + 0.0031(X> - 1.02). 



An alternative and more straightforward analysis of HT series is represented by the 
matching method, which was applied in Refs. 21,39 for the two- and three-dimensional Ising 
models. The idea is to generate sequences of estimates by fitting the expansion coefficients 
with their asymptotic form. By adding a sufficiently large number of terms one can make 
the convergence as fast as possible, although, of course, the procedure becomes unstable 
if the number of terms included is too large compared to the number of available terms. 
The estimates of the critical exponents from this analysis are consistent with the blAn 
results and with the fact that the 4 model at A ^ 2.1 and the ddXY model at D ?a 1.1 are 
approximately improved, i.e., the coefficients of the leading scaling corrections are consistent 
with zero. However, the results of this analysis are not sufficiently stable and precise to 
improve the estimates already obtained. 

Estimates of the critical exponents can be then obtained by evaluating Eqs. (88) and (89) 
at the FSS estimates of A*, i.e., A* = 2.15(5), and D*, i.e., D* = 1.06(2), where the residual 
effect of the leading scaling correction should vanish. We refer to these results as MC+IHT 
estimates. For the 4 theory we obtain 

v = 0.67170(4) [2]{9}, (90) 
7 = 1.31779(3)[5]{11}, 
r] = 0.03816(7){10} 

and using scaling and hyperscaling relations 

a = 2-3u = -0.01510(12)[6]{27}, (91) 
77 = 2 - 7/1/ = 0.03813(15){10}. 

The error due to the uncertainty on A* is reported in braces. For the ddXY model we obtain 

v = 0.67163(6)[2]{9}, (92) 
7 = 1.31774(11)[5]{14}, 
77 = 0.03811(20){9}, 
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and 



a = -0.01489(18)[6]{27}, 

,7 = 2-7/1/ = 0.03800(25){6}. 



(93) 



There is good agreement between the MC+IHT estimates obtained from the 4 and ddXY 
models. We stress that this represents a nontrivial check of the hypotheses underlying the 
IHT analysis. In particular, this supports our working hypothesis that effects due to next-to- 
leading nonanalytic corrections are negligible. Finally, estimates of the critical exponents 5 
and (3 can be obtained using the hyperscaling relations 5 = (5 — 77) /(H-77) and (3 = v{\+rj)/2. 

C. Combining IHT and FSS analyses 

Critical exponents can be also estimated by comparing Eqs. (88) and (89) with the analog 
FSS results (53) and (54) obtained from the MC data of Sjj 4 - We note that the coefficients 
that give the dependence on A and D in the IHT and FSS expressions have opposite sign. 
Therefore, the results of the two analyses agree only in a relatively small A (or D) interval. 
This is shown in Fig. 10. This comparison provides an estimate of v somewhat independent 
of the determination of the optimal values A*, D*. Taking as estimates of v values consistent 
with both the IHT and FSS analyses, we find the following FSS+IHT results 



The good agreement between the two results nicely support our estimates of the errors. The 
FSS+IHT results represent our most precise estimates of the critical exponent v. Corre- 
sponding estimates of a can be obtained using the scaling relation a = 2 — 3za 

Note that this analysis also provides alternative estimates of A* and D*. The intersection 
region in Fig. 10 indicates A* = 2.14(5) and D* = 1.07(3), which are in good agreement with 
those obtained in Sec. IV F. 

D. Universal amplitude ratios 

Using HT methods, it is possible to compute the first coefficients g 2 j and r 2 j appearing 
in the small-magnetization expansion of the Helmholtz free energy and of the equation of 



v = 0.67168(10) for the 4 model, 
v = 0.67167(15) for the ddXY model. 



(94) 
(95) 
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FIG. 10: Comparison of the HT and FSS results for the critical exponent v as functions of the 
parameters A, D around their optimal values. The dotted lines correspond to the estimates of v 
that are consistent with both the IHT and FSS analyses. 

state. 4 Indeed, these quantities can be expressed in terms of zero-momentum 2j-correlation 
functions and of the correlation length. They are denned as 

3 X4 

2 X 2 e 



34 = -o3rt> ( % ) 



and 



r 6 = 10-— —3-, (97) 
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etc... 

Using their extended HT series, we update the estimates of #4, r 6 , and r 8 obtained in 
Ref. 5. Consider a universal amplitude ratio A which, for t = 1 — (3/(3 c — ► 0, behaves as 

= A* + c^ A + c 2 t A2 + . . . + ai t + a 2 t 2 + ... (98) 

In order to determine A* from the HT series of A(t), we consider blAl's, whose behavior at 
(3 C is given by 

lAl*f((3)(l-(3/(3 c ) ( + g((3), (99) 
where f{j3) and g(/3) are regular at j3 c , except when ( is a nonnegative integer. In particular, 

In the case we are considering, ( is positive and therefore, g(f3 c ) provides an estimate of A*. 
Moreover, for improved Hamiltonians we expect ( = A 2 ~ 1.2, instead of ( = A ps 0.53. 

In the case of 54 we analyze the series /3 3//2 g 4 = X^i=o a */^- -'■ n Table VIII we report some 
results for several values of A and D. Assuming a linear dependence on A, D around their 
optimal values, we find 

g A = 21.153(4) + 0.14(A- 2.10) (101) 

for the 4 model, and 

#4 = 21.15(3) + 0A(D- 1.03) (102) 

for the ddXY model. We estimate the critical value of #4 by evaluating the above expressions 
at A* and D*. We obtain respectively 

g 4 = 21.160(4){7}, g 4 = 21.16(3){1}, (103) 

where the error in braces is related to the uncertainty on the estimates of A* and D*. We 
consider g± = 21.16(1) as our final estimate. This significantly improves our earlier result 5 
g 4 = 21.14(6) and is fully consistent with the FT estimate g± = 21.16(5) obtained from an 
analysis of six-loop perturbative series. 6 ' 40 Other results for g 4 can be found in Ref. 4. The 
results for the nonanalytic exponent (, reported in Table VIII, give 

C = 1.16(6) + 1.3(A -2.10) (104) 
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TABLE VIII: Estimates of the fixed-point value of g± from the 20th-order HT series of /3 3 / 2 <74. 



The error due to the uncertainty of (3 C is negligible. 

model <?4 C 

<p A , A = 2.07 21.149(5) 1.13(6) 

4 , A = 2.10 21.153(4) 1.16(6) 

4> 4 , A = 2.20 21.167(4) 1.28(7) 

ddXY,L> = 0.90 21.13(7) 0.9(2) 

ddXY, D = 1.02 21.15(3) 1.1(4) 

ddXY, D = 1.03 21.15(3) 1.1(3) 

ddXY, D = 1.20 21.22(2) 2.5(1.1) 



in the case of the 4 model. Evaluating ( at A = A* we obtain an estimate of A 2 , i.e., 
A 2 = 1.23(6){7}, corresponding to u 2 = 1.83(19). A consistent, but less precise, estimate 
can be obtained from the results for the ddXY model. 

It is worth mentioning some results obtained for the analysis of the HT series of g± in the 
standard XY model. Using blAl's, biased so that (3 C = 0.4541652(11) and A = 0.527(13) 
(corresponding to cu — 0.785(20)), we obtain g 4 = 21.12(5), which is in good agreement with 
the estimate obtained from the improved models, but much less precise. Moreover, if we 
analyze the same series biasing (3 C = 0.4541652(11) and g 4 = 21.16(1) and taking scaling- 
correction exponent A as a free parameter, we obtain the estimate A = 0.56(5), which is in 
agreement with the result obtained from the FSS analysis. 

Similar analyses applied to the 18th-order HT series of r§ and r§ provide the estimates 
r 6 = 1.96(2) and r 8 = 1.5(1), which substantially confirm those obtained in Ref. 5. These 
results can be used to compute approximations of the critical equation of state, using the 
method outlined in Refs. 5,41, which is based on an appropriate analytic continuation in 
the t, H space. Using our new estimates of a, r], r 6 and r 8 , we obtain results for the 
critical equation of state and universal amplitude ratios that are substantially equivalent 
to those obtained in Ref. 5, essentially because our new HT results do not significantly 
improve the estimates of r 6 , r 8 , and no precise and reliable estimates of the higher-order 
coefficients r 2 j are available. Therefore we do not provide further details. We only report 
the result R a = (1 — A + /A~)/a = 4.3(2) which is relevant for the superfluid transition 
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in 4 He. For comparison, we mention the FT result 42 R a = 4.43(8), the numerical MC 
results R a = 4.20(5) (Ref. 43) and R a = 4.0(1) (Ref. 44) and the experimental estimate 
R a = 4.154(22) reported in Ref. 3. 



1 J. A. Lipa, D.R. Swanson, J. A. Nissen, T.C.P. Chui, and U.E. Israelsson, Phys. Rev. Lett. 76, 
944 (1996). 

2 J.A. Lipa, D.R. Swanson, J.A. Nissen, Z.K. Geng, P.R. Williamson, D.A. Strieker, T.C.P. Chui, 
U.E. Israelsson, and M. Larson, Phys. Rev. Lett. 84, 4894 (2000). 

3 J.A. Lipa, J.A. Nissen, D.A. Strieker, D.R. Swanson, and T.C.P. Chui, Phys. Rev. B 68, 174518 
(2003). 

4 A. Pelissetto and E. Vicari, Phys. Rep. 368, 549 (2002). 

5 M. Campostrini, M. Hasenbusch, A. Pelissetto, P. Rossi, and E. Vicari, Phys. Rev. B 63, 214503 
(2001). 

6 R. Guida and J. Zinn-Justin, J. Phys. A 31, 8103 (1998). 

7 J. Zinn-Justin, Quantum Field Theory and Critical Phenomena, fourth edition (Oxford Univer- 
sity Press, Oxford, 2002). 

8 H. Kleinert and V. Schulte-Frohlinde, Critical Properties of cf) 4 Theories (World Scientific, Sin- 
gapore, 2001). 

9 M. Hasenbusch and T. T6r6k, J. Phys. A 32, 6361 (1999). 

P. Butera and M. Comi, Phys. Rev. B 60, 6749 (1999). 

1 M. Campostrini, A. Pelissetto, P. Rossi, and E. Vicari, Phys. Rev. B 61 5905 (2000). 

2 E. Burovski, J. Machta, N. Prokof'ev, and B. Svistunov, cond-mat/0507352. 

3 J.A. Lipa, S Wang, J.A. Nissen, and D. Avaloff, Space Advances in Space Research 36, 119 
(2005). 

4 A. Aharony and M.E. Fisher, Phys. Rev. B 27, 4394 (1983). 

5 J. Salas and A.D. Sokal, J. Stat. Phys. 98, 551 (2000); cond-mat/9904038vl. 

6 M. Caselle, M. Hasenbusch, A. Pelissetto, and E. Vicari, J. Phys. A 35, 4861 (2002). 

7 K.E. Newman and E.K. Riedel, Phys. Rev. B 30, 6615 (1984). 

8 H.G. Ballesteros, L.A. Fernandez, V. Martin-Mayor, and A. Munoz-Sudupe, Phys. Lett. B 441, 
330 (1998). 

43 



19 M. Hasenbusch, K. Pinn and S. Vinti, Phys. Rev. B 59, 11471 (1999). 

20 M. Hasenbusch, J. Phys. A 32, 4851 (1999). 

21 M. Campostrini, A. Pelissetto, P. Rossi, and E. Vicari, Phys. Rev. E 65, 066127 (2002); E 60, 
3526 (1999). 

22 M. Hasenbusch, J. Phys. A 34, 8221 (2001). 

23 M. Campostrini, M. Hasenbusch, A. Pelissetto, P. Rossi, and E. Vicari, Phys. Rev. B 65, 144520 
(2002). 

24 V. Privman, in Finite scaling and Numerical Simulations of Statistical Systems, edited by 
V. Privman, World Scientific, Singapore, 1990. 

25 M. Campostrini, A. Pelissetto, P. Rossi, and E. Vicari, Phys. Rev. E 57, 184 (1998). 

26 S. Caracciolo and A. Pelissetto, Phys. Rev. D 58, 105007 (1998). 

27 J. Salas, J. Phys. A 35, 1833 (2002). 

28 N.Sh. Izmailian and C.-H. Hu, Phys. Rev. E 65, 036103 (2002). 

29 U. Wolff, Phys. Rev. Lett. 62, 361 (1989). 

30 R.C. Brower and P. Tamayo, Phys. Rev. Lett. 62, 1087 (1989). 

31 M. Liischer, Comput. Phys. Commun. 79, 100 (1994). 

32 The G05CAF generator is a linear congruential generator: Ij + i = Mod[aIj + c, n] with a = 13 13 , 
c = and n = 2 59 . 

33 M. Creutz, Phys. Rev. Lett. 50, 1411 (1983). 

34 M. Creutz, Phys. Rev. Lett. 69, 1002 (1992). 

35 Y.J. Deng, H.W.J. Blote, M.P. Nightingale, Phys. Rev. E 72, 016128 (2005). 

36 M. Hasenbusch, A. Pelissetto, and E. Vicari, J. Stat. Mech. P12002 (2005). 

37 M. Campostrini, J. Stat. Phys. 103, 369 (2001). 

38 A.J. Guttmann, in Phase Transitions and Critical Phenomena, Vol. 13, edited by C. Domb and 
J. Lebowitz (Academic Press, New York, 1989). 

39 D. MacDonald, S. Joseph, D.L. Hunter, L.L. Moseley, N. Jan, and A.J. Guttmann, J. Phys. A 
33, 5973 (2000). 

40 In FT calculations, such as those reported in Ref. 6, one usually introduces the rescaled coupling 
g = 5<74/(247r). Our new estimate of 54 corresponds to g = 1.4032(7), to be compared with the 
FT result g = 1.403(3) (Ref. 6). 

41 M. Campostrini, A. Pelissetto, P. Rossi, and E. Vicari, Phys. Rev. B 62, 5843 (2000). 

44 



M. Strosser and V. Dohm, Phys. Rev. E 67, 056115 (2003). 

A. Cucchieri, J. Engels, S. Holtmann, T. Mendes, and T. Schulze, J. Phys. A 35, 6517 (2002). 
Preliminary result, M. Hasenbusch, in preparation. 



45 



